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failure properties, the burst pressure and fatigue lifetime. The favorable combination of high reliability 
and practicability of the composite vessel is a challenging task from the beginning of design. This paper 
gives a comprehensive review on recent development of numerical simulation and optimization for the 
designed composite vessel. First, methods on damage modeling for predicting the failure properties and 
degradation mechanisms of the composite vessel are reviewed. Second, research on predicting the burst 
pressure and lifetime of the composite vessel is reviewed. The academic work on the damage modeling, 
progressive failure analysis and finite element implementation which explains the failure properties and 
stiffness degradation mechanisms of the composite vessel is summarized. Computational methods on the 
burst pressure, the strength reliability and lifetime of the composite vessel are also evaluated. Finally, ideal 
design which aims to lessen the weight of a composite vessel to the maximum extent under strength 
and stiffness constraints is commented. The optimization efficiency using different algorithms is also 
comparatively studied. The numerical simulation and optimization as important fundamental research 
constitute a design platform for the composite vessel. It deserves pointing out the lightweight design 
conception as a remarkable tendency that combines advanced numerical methods and manufacturing 
technique develops rapidly, commits to improving the reliability and practicability of the composite 
vessel. It is expected the lightweight design technique plays an increasingly important role in developing 
the composite vessel as their value is further highlighted. 

© 2012 Elsevier Ltd. All rights reserved. 


Contents 
1;.. Introductions esii eenei pintre int a Yeas gees EAEE sa eet eaawelel ss gawd suadele Meade DEERE EED aba KAT EOE A Hebe dead conve ES EEE kE eiet 
2. Geometry structure of composite hydrogen storage vessel...... 


3. Finite element modeling of composite hydrogen storage vessel. 
Bul, -Design theories. ccsicesacs Gesesas rerun EEE ANNEER AOE 
3.1.1. Grid theory 

3.1.2. Composite laminate theory.....................008 

3.2. Finite element modeling 


4. Failure analysis and strength prediction of composite hydrogen storage Vessel ............cccccee cece cece cece cece eens ee enene eee eeeeeeeeeeeenenens 1821 
4:15 Damage constitutive modeling: rrsereoreserireres itensieden PENVEN EAEE EEE thes spain sd OEE EADE E CANET HENDE EEREN 1821 
4.2; Failüre criteria för Composites essenin rE rI E EED NENEA NEEE weed san REEE 302 deals E REO OASE bee 1821 
43: . Progressive failtire analysis escrocii erini rE EEE E EE EEE E E EEEE EAEE E E O TEE 1822 
AA... Burst pressure predictions mssiissis cosdcawane rinin Eras in EEEPC KEETE EEEE EENE RIE EEEE sees Reda mae HET EEIE CEE SEEE 1823 
5. Optimal design of composite hydrogen storage vessel...........sesesesesseesererseseseserererereeseseserorerseseseserererseseseeoeerereeseseeeerto 1824 
5i- Safety coefficient- based desig Mn 365 dike secctsns E r EEE TE dean RE E ate N E E N Fa ses 1824 
52: Reliability-and lifetime-based desioi ss sisic siiscccsines sieri ievos ro eaves E ET A TEE veld Vee E HEE ORE OE EN UTE TEE EES ETES 1824 


* Corresponding author. Tel.: +86 571 87953370; fax: +86 571 87953393. 
E-mail address: jyzh@zju.edu.cn (J.Y. Zheng). 


1364-0321/$ - see front matter © 2012 Elsevier Ltd. All rights reserved. 
doi:10.1016/j.rser.2012.01.006 


1818 P.F. Liu et al. / Renewable and Sustainable Energy Reviews 16 (2012) 1817-1827 


6. Further prospect of composite hydrogen storage vessel.... 
AcknowledgementS............cccseccesssesceettencceecoees 


Referentes 525 sasiansad dane natbe bined vee cadulenemanbee dee pabaaees sae dead eek ce 


1. Introduction 


The energy crisis and environment pollution all around the 
world enable it necessary to explore advanced energy. Hydrogen 
energy is commonly considered as important renewable and sus- 
tainable energy in the new century because of a series of advantages 
such as cleanness, high-efficiency and fruitful source [1-4]. Hydro- 
gen can be generated by solar energy, wind energy, biomass and 
water [5-7]. Yet, the use of hydrogen energy deals with many tech- 
nique problems such as the hydrogen generation, the hydrogen 
storage and transport, and the hydrogen combustion. To promote 
the hydrogen energy industry and potential hydrogen economy, it 
is critical to develop hydrogen energy technique [8-10]. 

The hydrogen fuel cell car is typically representative of devel- 
oping hydrogen industry. Safe, high-efficiency and economical 
hydrogen storage technique is a key to ensure favorable run of 
hydrogen energy car. Among many hydrogen storage patterns 
including high-pressure gaseous storage, cryogenic liquid stor- 
age and chemical hydrogen storage, high-pressure gaseous storage 
has become the most popular technique [11-14]. For example, 
the hydrogen fuel cell cars used in both 2008 Beijing Olympic 
Games and 2010 Shanghai World Exposition used composite high- 
pressure storage vessel. 

Because of the limited space available for the car, the basic 
design requirements for the storage vessel are safety, reliability 
and economy. However, the composite pressure vessel may work 
under the high-pressure and high-temperature environment. This 
not only presents a strong challenge to the physical and mechan- 
ical performance, but also to the reliable and economical design 
about how to achieve a perfect combination of safety performance 
and low cost [15]. Conventional movable metal pressure vessel can 
no longer be competent for the rigorous need of high strength - 
and stiffness - weight ratios. In this case, the composite filament 
wound technology is introduced to improve the storage efficiency 
[16-18]. To reach the aim of large-capacity hydrogen storage, the 
composite vessel includes an inner aluminum layer and se outer 
carbon fiber/epoxy composite wound layers in China [19]. In gen- 
eral, the composites are used to make the pressure vessel by placing 
them in different orientations for different layers and ina common 
orientation within a layer. These layers are stacked in such a way 
to achieve high stiffness and strength [20]. The designable proper- 
ties of the composite vessel enable it flexibly applicable to different 
working conditions and environments by adapting the fiber wound 
patterns, the dome shapes and the layer thickness. 

The design of the composite vessel as a fundamental research 
work relates the physical and mechanical properties of materials 
to the geometry. Essentially, two main works are required for the 
design: one is the prediction of burst pressure and another is the 
optimization. The burst pressure denotes the limit load-carrying 
ability of the composite vessel. Since the failure mechanisms of 
composites are complicated from the view of composite microme- 
chanics, including the matrix cracking, fiber/matrix debonding and 
fiber breakage as well as their interactions [21,22], accurate pre- 
diction of the burst pressure is a challenging work. Besides, proper 
explanation of the stiffness degradation as the physical response to 
the continuous failure also needs special theory. 

With the development of computer, the finite element analysis 
has become a powerful tool to deal with the numerical problems 
of complex structures [23]. By combining the continuum damage 
mechanics (CDM) and the finite element analysis (FEA), the pro- 
gressive failure analysis of the composite vessel is used to account 


for the stiffness degradation and to predict the burst pressure 
[24,25]. Especially, the FE technique that implements the damage 
model arouses wide academic attention. 

In terms of the optimization of composite vessel, several intelli- 
gent algorithms that simulate the evolutionary principles in the 
disciplines of computational biology, physics and immunology, 
have been used. The popular algorithms are the genetic algorithm, 
simulated annealing and artificial immune system [26-28]. The 
optimization attempts to improve the weight, strength, reliability 
and lifetime of the composite vessel by designing the wound angles 
and thickness. 

As the design method and manufacturing technique advance, 
the composite vessel takes on the approach to the lightweight 
development. The weight decrease has become a technique bot- 
tleneck which restricts the economy and practicability of the 
composite vessel. In this case, it is urgent that a set of material- 
structure-processing integrated design and calculation method of 
the composite vessel be established. It is true there are some other 
manufacturing and test reasons affecting the performance of com- 
posite vessel, so what a large role the design plays may require 
evaluated. However, the importance of numerical simulation and 
optimization is expected to be further highlighted by the present 
review. 


2. Geometry structure of composite hydrogen storage 
vessel 


Fig. 1 shows a composite vessel used in the hydrogen fuel car. 
Fig. 2 shows the basic structure and wound pattern of composite 
vessel. The vessel includes an inner aluminum layer and several 
outer composite layers. The aluminum layer prevents the gas leak- 
age and provides a mould for filament wound, and the composite 
layers are responsible for resisting the internal pressure. For the 
cylinder part, spiral wound and hoop wound patterns are used. Yet, 
only spiral wound pattern appears at the head part. 

From the manufacturing principle, the wound trace for compos- 
ite layers is a combination of the rotation of aluminium liner and 
the axial movement of wound machine. The spiral wound angle ap 
at the cylinder is determined by [29] 


a = arcsin (7) (1) 


0 


Fig. 1. Lightweight composite high-pressure hydrogen storage vessel used in the 
hydrogen fuel-cell vehicle. 
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Fig. 2. Basic structure and wound pattern of composite vessel: (a) hoop wound and 
(b) sprial wound. 


where r is the radius of the polar axis and Rọ denotes the inner 
radius of the cylinder. 

The spiral wound at the head obeys the geodesic path algorithm, 
marking the shortest distance between any two points on the head. 
The spiral wound angle « at the head changes from 90° at the polar 
axis to @p at the cylinder [29] 


a = arcsin (7) (2) 


The thickness H of the composite layer at any radius R of the 
head is calculated by [29] 


(R2 — r2) 
(R= 7?) 
where h is the wound thickness at the cylinder. From Eqs. (1)-(3), 


the polar radius r affects the wound angle «œ and the thickness H at 
the head. 


I 
ll 
> 


(3) 


3. Finite element modeling of composite hydrogen storage 
vessel 


Accurate and fast modeling of the composite vessel is an impor- 
tant work in the design, which involves many parameters such as 
the wound angle and thickness. Parametric modeling by exacting 
these parameters as the design variables facilitates the strength 
prediction and optimization for the composite vessels with differ- 
ent sizes. 


3.1. Design theories 


In general, the composite vessel is considered as the composite 
laminated structure. Now, there are two main design theories: grid 
theory and composite laminate theory. 


3.1.1. Grid theory 
The basic assumptions according to the grid theory are: (1) only 
the longitudinal carbon fiber bears the pressure and (2) effects of 


wound patterns are neglected. The grid theory can be used to cal- 
culate the thickness of each composite layer and the longitudinal 
in situ fiber strength. The strength value is inversely solved by the 
burst pressure of composite vessel. Chen [30] gave a set of design 
method of the composite vessel based on the grid theory, and the 
total composite wound thickness hy is expressed as 


hf = Roh ; hø = RoPb o2 — tan? ag) 
20) COS? ag 20% (4) 
h Haah 3RoPp 
= + hø = 
p= Nya + he = See 


where hy, and hy are longitudinal and hoop wound thickness, Pp is 
the burst pressure of composite vessel, and opp is the in situ strength 
of carbon fiber. 

It should be pointed out the grid theory is merely an ideal design 
method which largely depends on the processing parameters and 
the sizes of the experimental composite vessel. This shows seeking 
a more efficient design theory is necessary. 


3.1.2. Composite laminate theory 

The basic assumptions of the classical composite laminate the- 
ory (CCLT) are: (1) perfect bonding appears at the interface between 
each layer, (2) the mechanical properties of composite laminates 
are substituted by those of a middle plane, and (3) the normal 
stress is neglected at the section parallel to the middle plane. Lif- 
shitz et al. [31] proposed a method using the CCLT to calculate the 
stress and strain in non-symmetric wound pressure vessel with 
thick metal liners. However, the last two assumptions above do 
not strictly hold for 3D composite vessel. In this case, Parnas et al. 
[32] derived the elastic stress and displacement solutions for 3D 
cylindrical composite laminates based on the CCLT and general- 
ized plane strain assumption. As the liner is assumed isotropic 
and the composite layers are considered transversely isotropic, 
Chapelle et al. [33] further derived the analytical solutions of 
the stress, strain and displacement for 3D cylindrical composite 
laminates. 

Fig. 3(a) shows the cylindrical part of composite vessel and 
Fig. 3(b) describes a representative volume element taken from 


r 2 
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liner 
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Fig. 3. (a) Composite hydrogen storage vessel (transverse section) and (b) schematic 
illustration of the relationship between the on-axis coordinate (1,2,3) and the off- 
axis coordinate (r,9,z). 
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Fig. 3(a), showing the principal direction (1,2,3) of acomposite layer 
under the cylindrical coordinate (r,9,z). If the axial strains at all lay- 
ers are assumed constant and the shear strains are independent of 
z, the strain-displacement relationship is given by [33] 


k k k (k) 
do dP w dy ae tay M 
r dr’ 6 r’ Zz dz 0; 
(5) 


where ez is the axial strain, yz is the shear strain and yo is twist 
per unit length. 

The off-axis stress-strain relationship under the hygro-thermo- 
mechanics multiphysics field is expressed as 


(k) = = a a (k) (k) (k) (k) 


where the symbols Õli, j = 1, 2, 3, 6) are off-axis elastic con- 
stants of composites. The symbols oz, og, or and Tz are the off-axis 
axial, hoop, radial and shear stresses, respectively. £z, £9, €r and yz9 
are the corresponding strains, respectively. AT and AC are tem- 
perature and moisture concentrations, respectively. œ;(i= 1, 2, 3, 6) 
and §;(i=1, 2, 3, 6) are the corresponding thermal and moisture 
expansion coefficients. 

For the kth layer, the equilibrium equation under the cylindrical 
coordinate is given by [33] 


1arot) oH" _ (7) 
ro r 


After introducing the displacement continuity conditions, the 
radial stress continuity conditions as well as the axial equilibrium 
and zero torsion conditions between two neighbouring layers, the 
analytical solutions for the strain, stress and displacement can be 
derived by Eqs. (5)-(7). 


3.2. Finite element modeling 


The finite element analysis and CCLT are combined to consti- 
tute the design platform since the analytical solutions cannot be 


competent for progressive failure analysis. Yan et al. [34] discussed 
the finite element technique of the composite vessel including finite 
element modeling and commercial finite element codes. Besides, 
the reliability, efficiency and cost as well as the further prospect 
are discussed. Liu and Zheng [35,36] established parametric finite 
element model for the cylinder part using ANSYS finite element 
code. Wang et al. [37] established the finite element model of 
the composite vessel with variable thickness and angles at the 
head. Li et al. [38] established the design model of the compos- 
ite vessel using ANSYS-APDL (Ansys Parametric Design Language) 
and VC++ language. However, the design platform has not been 
completely established since ANSYS software does not provide an 
efficient modeling module for composites until now. Based on the 
composite micromechanics and modified CCLT, Antunes et al. [39] 
established the finite element model of the composite vessel using 
ABAQUS and Algor finite element codes. Park et al. [40] calculated 
the filament wound patterns using semi-geodesic path equations 
for an arbitrary surface, and established the finite element model of 
the composite vessel using ABAQUS finite element code. Further, 
Kim et al. [41] used the semi-geodesic path algorithm to calcu- 
late possible wound patterns by considering the windability and 
slippage between the fiber and the mandrel surface. Yet, these 
work concentrates on the specific composite vessel and does not 
perform parametric modeling. Recently, we established the para- 
metric finite element model by combining the ANSYS and ABAQUS 
finite element codes. Fig. 4 shows the flow chart of finite ele- 
ment modeling of composite vessel, and Fig. 5 shows a 90° mesh 
model. A subroutine is written using ANSYS-APDL to establish the 
geometry and mesh model, where some characteristic parame- 
ters such as the composite lay-up thickness, the cylinder inner 
radius and thickness, the liner thickness as well as the mouse 
radius and thickness are extracted as the design variables. Then, 
the mesh model is imported into ABAQUS by the ANSYS/ABAQUS 
file transformation interface. After that, a local cylindrical coor- 
dinate system (r,0,z) is defined in ABAQUS and the composite 
lay-up module is employed to assign the lay-up angles and mate- 
rial properties for composite elements. Fig. 6 shows schematic 
stress distributions under working pressure for a composite ves- 
sel. The wound angles at the cylinder are 90°, —90°, ao, —ao, 90°, 
—90°, œo, —&o, 90° and —90° in turn. It can be seen that the maxi- 
mum principal stress appears at the joint of the head and cylinder 


Fig. 4. Flow chart of finite element modeling of composite vessel. 
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Fig. 5. Parametric finite element model of composite vessel. 


for the spiral wound layer, in good agreement with Yang’s results 
[29]. 


4. Failure analysis and strength prediction of composite 
hydrogen storage vessel 


Currently, finite element technique has been widely used to 
predict the failure properties and strengths of composites. In the 
following, some recent research on the failure analysis and strength 
prediction of the composite vessel is reviewed. 


4.1. Damage constitutive modeling 


With the accelerating application to the hydrogen fuel cell car in 
China, the composite vessel may be subjected to the hygro-thermo- 
mechanics environment arising from the subtropical climate. The 
aging of epoxy matrix leads to the stiffness and strength degrada- 
tion, thereby affecting the lifetime of composite vessel. In this case, 
it is significant to explore the failure properties using progressive 
failure analysis. 

Strictly speaking, there are two distinct failure mechanisms 
for the composite vessel: intralaminar damage and interlaminar 
delamination. From the burst experiments of composite vessels, 
the intralaminar damage instead of interlaminar delamination 
affects the mechanical performance of the composite vessel. This 
is because two neighbouring composite layers still keep normally 
pressing and tangentially slipping despite a little interlaminar fail- 
ure appearing before the burst [42]. Thus, only the intralaminar 
failure is generally considered. 

Since Kachanov [43] proposed the CDM to study the creep rup- 
ture of metals, progressive failure analysis using the CDM has 
been widely used to predict the stiffness degradation and damage 
evolution properties of composites by defining phenomenological 
damage tensors and damage evolution laws [44-50]. A schematic 
3D damage model from Eq. (6) under the hygro-thermo-mechanics 
field is expressed as 


eal hie éz—a, AT- f Ac] (—dp)Cn 
og -c| £0- %2AT — Bz AC c, = | = dA = dm)C21 
Or =d] 6, —a3AT — B3 AC ae (1 — dp )Ci3 

T29 Yz0 — AG AT — Bg AC Ge 


where d,(i=f, m, s) are the damage variables representing the 
fiber breakage, matrix cracking and fiber/matrix interface failure, 
respectively. 


According to Lemaitre and Chaboche [51,52], Ristinmaa and 
Ottosen [53], Schapery [54], Murakami and Kamiya [55], Hayakawa 
et al. [56] and Simo et al. [57], the damage constitutive modeling 
based on the CDM is summarized. For the damaged materials, the 
Cauchy stress tensor o can be substituted by the nominal stress 
tensor o 


o = (1 — D)ō (9) 


where 0<D<1. D=0 represents the perfect materials and D=1 
denotes the completely damaged materials. 

The Helmholtz free energy per unit mass W for elastic materials 
under isothermal conditions is written as 


pw = pw'(eq— ef, Dy) + pW) (10) 


where w® and y are the free energy which represents the elastic 
deformations and damage hardening, respectively. « is an internal 
variable and p is the density. 

The thermodynamic conjugate forces (Y;;, B) corresponding to 
the internal variables (Dj, «) are expressed as 


x a1) 


The Clausius-Duhem dissipation inequality according to the 
CDM is written as 


Y = -pý + oye; > 0 (12) 
where y is the power of dissipation due to damage. If the dissipation 
y reaches the maximum, the damage evolution law is given by 

. q OF4 jd oF@ 
4 OY,’ ðB 


(13) 


After introducing three damage potential functions FUY;, D, B) 
for three intralaminar failure mechanisms for composites, the con- 
sistency parameter àf, the damage variable D and elastic matrix C4 
are obtained as the functions of stress and strain. 


4.2. Failure criteria for composites 


Failure shows the degradated ability to resist the external 
load. The failure properties of a composite vessel depend on the 
geometry sizes, the material features, the loading history and envi- 
ronment. Over the past five decades, the popular failure criteria 
are the maximum stress, Hashin, Hoffman, Hinton, Christensen, 
Rotem, McCartney, Yamada-Sun, Chang-Chang, Tsai-Hill, Tsai-Wu, 
Puck and Huang’s criteria [58-70]. In terms of the tensor notation, 
the generalized Tsai-Wu criterion is written as 


Fio; + Fyojoj + Fijkoi0joęk = 1, i,j,k=1,2,...,6 (14) 


where F is strength coefficient tensor. 

Such as the Tsai-Wu, Tsai-Hill and Hoffman failure criteria 
identify the failure of composites, but they cannot identify the fail- 
ure modes. In contrast, the maximum stress, Hashin, Christensen, 
Rotem, McCartney, Yamada-Sun, Hinton, Chang-Chang, Puck and 
Huang criteria can identify the failure modes by defining the fail- 
ure criteria for the tensile and compressive failure as well as the 
fiber breakage, matrix cracking and delamination. 


-d1 -dm O-d)Č3 Ge = 
(1 — dm)C22 3 C26 (8) 
C23 C33 C36 | 
C26 C36 (1 — ds )dsCe6 


Hinton et al. [60] explored existing failure theories for com- 
posites by launching a world-wide failure exercise, and tested 
some leading failure theories by comparing experimental data. 


1822 P.F. Liu et al. / Renewable and Sustainable Energy Reviews 16 (2012) 1817-1827 


#116 7e+02 
*1.006e+02 


(c) 


S, Max. Principal 


Fig. 6. Schematic stress distributions at (a) liner, (b) hoop wound layer and (c) spiral wound layer for a designed composite vessel. 


Nahas [71], Tsai [72] and Sleight [73] performed extensive liter- 
ature surveys on the failure theories of composites and classified 
the failure criteria into two types: the non-interactive failure 
criteria and interactive failure criteria. The non-interactive fail- 
ure criterion is defined as one without interactions between the 
stress or strain components. These criteria compare the individ- 
ual stress or strain components with the corresponding material 
allowable strengths, and indicate the type of failure modes. The 
maximum stress, Hashin, Christensen, Rotem, McCartney, Yamada- 
Sun, Chang-Chang, Puck and Huang’s criteria belong to this type. 
In contrast, the interactive failure criteria consider interactions 
between stress components. These theories use a polynomial based 
on the material strengths to describe a failure surface. The Tsai-Wu, 
Tsai-Hill and Hoffman failure criteria belong to this type. Currently, 
there is also an evolving trend to develop the fracture energy-based 
failure criteria to predict the failure of composites using such as 
the cohesive fracture energy-based criterion [74,75] and finite ele- 
ment failure-based criterion [76,77]. However, accurate prediction 
of fracture properties of composites requires advanced experimen- 
tal technique. 


4.3. Progressive failure analysis 


Essentially, progressive failure analysis of the composite ves- 
sel is to implement the proposed damage constitutive model, 
the failure criteria and damage evolution law. Garnich and Akula 
[78], Orifici et al. [79], and Liu and Zheng [24] gave reviews 
on the current methods on the progressive failure analysis and 
finite element implementation of composites. The finite ele- 
ment equation under the hygro-thermo-mechanics field is given 
by 


Ku = Fs + Fr + Fc (15) 


where K=f vB'C,B dV is the stiffness matrix of a finite element, 
Fs= f sN'F; dS, Fr= f yB'Cgo ATdVand Fc= f yB'CyB ACdVare ele- 
ment forces arising from the surface load, thermal expansion and 
moisture expansion, respectively. B=LN! is strain matrix, L is dif- 
ferential operator and N is shape function. 

The flow chart for the progressive failure analysis of the com- 
posite vessel using finite element method is illustrated in Fig. 7, as 
summarized by five points [73]: (1) for each load step, the finite 
element analysis is performed and the on-axis stresses/strains at 
each element are obtained, (2) the stresses/strains at each ele- 
ment are compared with the material allowable values and used 
to determine whether some elements fail according to the failure 
criteria. If no failure is detected, the applied load is increased and 
the analysis continues, (3) if some elements fail, the damage vari- 
ables are updated, and the element stiffness constants are degraded 
according to the damage model. After that, the equilibrium of the 
composite vessel is re-established using the modified stiffness. This 
adjustment accounts for the nonlinear solution due to the stiffness 
change, (4) the calculations are performed continuously under the 
same load until there is no failed element anymore, (5) the iterative 
process repeats until the catastrophic failure showing the burst of 
the composite vessel. 

Chang [80] and Ju et al. [81] performed the initial and progres- 
sive failure analysis of the composite vessel based on the CCLT. 
However, the analysis is only limited to the analytical solution. 
Further, Perreux and Thiebaud [82], Liu and Zheng [35] per- 
formed the progressive failure analysis of the composite vessel 
by proposing damage models and developed the corresponding 
finite element technique. The numerical convergence problem is 
solved by introducing the arc-length algorithm or viscous sta- 
bilization method [83,84]. The fracture localization problem is 
addressed by adding a characteristic parameter into the damage 
model to eliminate the mesh effect [85]. However, to develop the 
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Fig. 7. Flow chart of progressive failure analysis of composite vessel. 


finite element technique for the failure analysis under the hygro- 
thermo-mechanics multiphysics field may still require a long 
way. 

Since the failure of composite is the macroscopic representation 
of the microscopic fiber/matrix failure evolvement, it is insufficient 
to explore the failure properties of composite vessels merely by 
using the CDM theory. In the past two decades, there were continu- 
ous efforts for developing hierarchical and multiscale finite element 
methods to perform the progressive failure analysis of composites. 
The multiscale damage model at the microscopic and macroscopic 
scales for the composite vessel is shown in Fig. 7. Square and hexag- 
onal cell models are established at the microscopic scale based on 
the observation of carbon fiber/epoxy composites using scanning 
electron microscope. Thus, the macroscopic damage tensors above 
can be represented by microscopic damage information. He and 
Wang [86] studied the transverse matrix cracking and stress con- 
centration for the composite vessel by multiscale failure analysis. 
Haj-Ali and Muliana [87] proposed an integrated micromechan- 
ics framework for the nonlinear analysis of composites. However, 
these work is only limited to the simple composite structures. Fur- 
ther, Blassiau et al. [88] established the microscopic cell models 
on the fiber/matrix load transfer and the multiscale model on the 
progressive failure analysis of the composite vessel. However, the 
failure properties of the microscopic epoxy matrix was not con- 
sidered and the damage initiation information for the fiber/matrix 
interface cannot be acquired. With the development of computer, 
the computational ability for the multiscale analysis linking the 
microscopic to macroscopic failure analysis is largely improved. 
Potential difficulty lies how to reveal all microscopic failure mech- 
anisms ina realistic way which may be elusive of engineering tests 
[74]. 


4.4. Burst pressure prediction 


Burst pressure of the composite vessel is the ultimate aim of 
progressive failure analysis. Mao et al. [89] estimated the burst 
pressure based on the statistical approach and experiments. Cohen 
et al. [90,91] investigated the effects of manufacturing and design 
variables such as the material parameters and wound parameters 
on the strength and stiffness of the composite vessel by exper- 
iments. Hwang et al. [92,93] performed the probabilistic failure 
analysis and experimental tests to predict the size effect on the 
burst pressure of composite vessels using Monte Carlo method. 
However, no information on the progressive failure properties of 
composites is obtained from these work. Sun et al. [94] performed 
the nonlinear finite element analysis to predict the burst pressure 
of the composite vessel using the maximum stress failure criterion 
and a stiffness degradation model. Onder et al. [95] explored the 
influence of temperature and wound angle on the burst pressure of 
the composite vessel by the modified CCLT and experiments. How- 
ever, they defined only the constant degradated stiffness regardless 
of the evolvement of damage variables. Based on the modified CCLT 
and CDM, Liu and Zheng [35,96] proposed damage models and finite 
element algorithm, respectively, to predict the burst pressure of the 
composite vessel by considering the progressive failure properties 
and nonlinear damage evolvement. Prediction of the burst pres- 
sure is essentially related to the numerical convergence problem 
in the finite element analysis. According to the design requirement 
of the composite vessel, only little matrix cracking under work- 
ing pressure is permitted. As the pressure increases, progressive 
interface debonding and matrix cracking appear. Only when a large 
number of fiber breakage suddenly appears is the composite vessel 
considered to burst. 
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5. Optimal design of composite hydrogen storage vessel 


The aim of design for the composite vessel is to reach the 
best match between the reliability and practicability. Advanced 
design technique can shorten the duration and save the cost. Now, 
the ISO (International Standardization Organization) and European 
Union have issued the standards for the design of composite ves- 
sel [97,98]. Unluckily, there is still no mature national standard of 
composite vessel until now in China since most of the design and 
manufacturing depend on the experience and experiment [99,19]. 
In this case, developing a set of efficient design method is nec- 
essary. Now, two main design methods for the composite vessel 
sprout: safety coefficient-based design, and reliability and lifetime- 
based design, which are considered as the improvement over the 
conventional experiential design. 


5.1. Safety coefficient-based design 


According to Leung [100], the working pressure of composite 
vessel is taken as the burst pressure divided by the safety coeffi- 
cient. However, the difficulty lies in how to select a proper safety 
coefficient which yet depends on the design and manufacturing 
levels. The CGH2R Draft Revision 10 of European Union and the 
ISO/DIS 15869.2 specify the safety coefficients of 2.35 and 2.25, 
respectively [97,98]. Then, the optimal design attempts to reach 
the minimum weight and highest strength under the constraint of 
safety coefficient. Tummala et al. [101] performed the multiobject 
optimization of the composite vessel on the pressure, temperature 
and weight. Krikanov [102] proposed a new optimization method 
to design the composite vessel with the minimum mass under 
strain and strength constraints. Especially, Liang et al. [103] car- 
ried out the optimal design of dome contours for the composite 
vessel subjected to geometric limits and wound condition. Wang 
and Zheng [37] proposed a mathematical design model based on 
the finite element analysis. However, the optimization efficiency 
above is low since only simple optimization methods are used. In 
this case, searching advanced intellectual algorithms becomes nec- 
essary. Ghiasi et al. [104] compared the optimization efficiency of 
various methods such as the gradient-based method, direct search 
method, genetic algorithm and simulated annealing. Messager et al. 
[105] developed a genetic algorithm-based optimization method 
to maximize the stability buckling pressure of the composite ves- 
sel. Kim et al. [106] proposed an optimal design method of the 
composite vessel by associating the finite element analysis with 
the genetic algorithm. Liu et al. [107,108] proposed an adaptive 
genetic algorithm and an immune algorithm, respectively, to per- 
form the optimal design of the composite vessel which aims to 
reach the minimum weight under the constraints of safety coef- 
ficient and processing parameters. The parametric model for the 
composite vessel is shown in Fig. 8. The intellectual optimization 
algorithms show higher convergence speed and precision than the 
gradient-based method. 


5.2. Reliability and lifetime-based design 


For a designed composite vessel, economy and practicability 
are based on safety and reliability. The reliability and lifetime of 
the composite vessel relate to many factors such as the material 
parameters, wound patterns and manufacturing defects. It is not 
always feasible to expect to use a single safety coefficient to pro- 
vide the guarantee for any occasional accident. Thus, it is necessary 
to develop the reliability and lifetime-based methods to perform 
the probabilistic design. Compared with the deterministic design 
method, the probabilistic design method may avoid costly over- 
design and conservative manufacturing. 
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Fig. 8. Optimization model of composite vessel. 


Now, the reliability and lifetime-based design conception has 
been successfully applied to the composite vessel. Béakou et al. 
[109,110] optimized the wound angle of cylindrical composite 
laminates using reliability methods for the composite structure. 
Gomes et al. [111] performed the optimization of laminated com- 
posites under the reliability constraint using genetic algorithm and 
neural network. However, the geometry model in there work is 
much simpler than the composite vessel. Further, Shen et al. [112] 
launched the reliability-based probabilistic design of the compos- 
ite vessel by combining the statistical theory and experiments. 
However, the probabilistic design is only limited to the analytical 
solution using the grid theory. Cai et al. [113] developed a set 
of reliability-based design method of the composite vessel and 
performed the sensitivity analysis on the influence of statistical 
variables on the safety factors under external pressure. Based on 
the progressive failure analysis, Liu and Zheng [114] performed 
the strength reliability-based design of the composite vessel using 
Monte Carlo simulation and response surface method, in which 
the lay-up thickness and angles of composite layers are consid- 
ered as the stochastic variables. Besides, Wang [115] compared 
several existing lifetime prediction technique of the composite 
vessel under static and cyclic internal pressure. Camara et al. 
[116] studied the strength reliability, the failure probabilities and 
the lifetime of the composite vessel by multiscale modeling and 
composite micromechanics. Bie et al. [117] proposed a theoretical 
model to predict the fatigue lifetime of the composite vessel by 
combining the micromechanics and CDM. However, the research 
above did not touch upon the lifetime-based design. Comond et al. 
[118] proposed a set of fatigue lifetime-based design method of 
the composite vessel based on the modified CCLT and the CDM. 
However, the design for the head part was not considered in 
their work. Further, Bertin et al. [119] explored the effects of 
wound patterns on the fatigue lifetime of the composite vessel 
under thermomechanical cyclic loading by experiments. It can be 
concluded the strength reliability and lifetime-based material- 
structure-processing integrated design technique will become an 
important tendency for the design of composite vessel. The best 
design method is to set up the high-fidelity numerical simulation 
and advanced optimization platform. But, the design must also 
consider too many factors including the statistical properties of 
fiber strengths, complex geometry structure, and efficiency and 
stability of numerical calculations. For the probabilistic design, the 
prediction of the burst pressure and fatigue lifetime of composite 
vessel is main limitation which otherwise needs more funda- 
mental research. After all, the complicated failure mechanisms of 
composites are still beyond our understanding until now [120]. 
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6. Further prospect of composite hydrogen storage vessel 


In recent years, the composite pressure vessels develop rapidly 
in the field of hydrogen storage. Academic and engineering atten- 
tions are paid to the design and manufacturing technique of them. 
With the development of computer-aided design technique, the 
composite vessel offers designers much design space. The design is 
to reach the ideal match among the reliability and lifetime, weight 
and cost. 

It is commonly recognized the composite vessel should take on 
the lightweight road [121]. The lightweight design, manufacturing, 
test and evaluation technique guarantee the safe and economical 
application of composite vessel. The ideal design method is to set up 
a set of material-structure-processing integrated design platform, 
and to solve the numerical simulation and optimization problems 
efficiently. This is not easily realized in view of the complicated 
relationships between the design requirements and design param- 
eters. 

The composite vessel is an important application of such as 
the advanced carbon fiber/epoxy composites to the area of hydro- 
gen energy. The complicated failure mechanisms and degradation 
mechanisms are distinct characteristic of composites although they 
exhibit high stiffness- and strength-density ratios. Further, the 
hygro-thermo-mechanics environment is also a large problem in 
the design of the composite vessel. More fundamental research on 
the failure properties of composites under the multiphysics field 
remains to be performed. 
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